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Abstract. We present VLBI and MERLIN data at 1.66 
and 4.99 GHz on the central component coincident with 
the nucleus of the giant radio galaxy, 3C236. The nuclear 
radio structure is composed of two complexes of emission 
which are resolved on scales from 1 milli-arcsec (mas) to 
1 arcsec. Oscillations with an amplitude of ~ 5° can be 
seen in the compact radio structure. Spectral index dis- 
tributions are plotted at angular resolutions of 10 and 25 
mas and allow us to identify the core component in the 
south-east emission complex. Re-examination of the HST 
WFPC-2 image of 3C236 by de Koff et al. ~ 



200C 



shows 

that the normal to the dust disk in the nucleus is ~ 30° 
from the plane of the sky and within 12° of parallel to the 
overall orientation of the radio source. We suggest that 
the radio axis is also at an angle of ~ 30° to the plane of 
the sky and that the north-west jet is on the approaching 
side. This orientation implies an overall size of 4.5 Mpc 
(H = 75 kms^Mpc" 1 , q Q = 0.5) for 3C236. The coin- 
cidence of a dust feature and the south-east compact jet, 
within the astrometric errors, leads us to suggest that the 
dust may be in the form of a cloud encountered by the 
jet in the first ~ 400 pc of its journey out from the nu- 
cleus. One-sided emission at 5 GHz on 1 mas scales would 
suggest that the jets are ejected initially at < 35° to the 
line of sight, but this is difficult to reconcile with the ob- 
vious orientation stability of the jet system as a whole. 
Free-free absorption of the counter-jet may be an alter- 
native explanation for the one-sideness. At the resolution 
of WSRT data at 327 MHz, the jet to the south-east is 
apparently continuous over a distance of 2.5 Mpc, mak- 
ing this the largest jet known in the universe. It is likely, 
however, that activity in the nucleus of 3C236 is episodic 



but with a shorter duty cycle than in the double-double 
source s stud ied by Schoenmakers et al. ( 2000 ) and Kaiser 
et al. d2000| ). 



Key words: Galaxies: nuclei - galaxies: jets - galax- 
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1. Introduction 

Extragalactic radio sources have been studied for many 
years, but it is still unclear how they are formed and how 
they evolve. A crucial element in the study of their evolu- 
tion is to identify the young counterparts of 'old' Fanaroff- 
Riley (FR) clas s 1 an d 2 extended objec ts. Re cent analyses 
by Fa nti et al. ( |1995| ), Readhead et al ( |1996[ ) , and Snellen 
et al. ( 200C| ) suggest that good candidates for young radio 
sources are the Gigahertz Peaked Spectrum (GPS) and 
Compact Steep Spectrum (CSS) sources with symmetric 
morphologies. They are small in angular size as expected 
for young sources, and could plausibly expand to larger 
sizes while decreasing in luminosity as required to match 
the statistics of the numbers of compact and extended 
objects. The propagation velocities of hot-spots in some 
small size symmetric sources have been measured at a few 
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tenths of the speed of light by Owsianik & Conway (1998) 
and Tschager et al. (200C); these velocities indicate dy- 
namic ages of a few hundred to a few thousand years. In a 
few objects we see both a GPS or CSS component associ- 
ated with the nucleus of the galaxy and extended FR1 or 
FR2 lobes either side of the nucleus. These sources with 
recurrent activity provide us with case studies of the evo- 
lution of individual objects which constrain models for the 
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evolution of radio sources in general. One such object is 
3C236. 

3C236 is the largest radio source known in the uni- 
verse and this in itself makes the source a powerful lab- 
oratory for the study of radio source evolution. The re- 
markable agreement between the orientation of the inner 
and outer radio structure was one of the original argu- 
ments that a massive black hole was the source of the 
AGN energy. Its angular extent of 39 arcmin corresponds 
to a projected size of 3.9 Mpc at the redshift of its asso- 
ciated 17 m galaxy (z = 0.1005 ± 0.0005, Hill et al. |1996 , 
H = 75 kms" 1 Mpc" 1 , q Q — 0.5). Its total radio luminos- 
ity of 1.5 x 10 43 ergs" 1 (Mack et al. |1998|) puts it in the 



transition zone between FR1 and FR2 systems. Its mor- 
phology is also intermediate between FR1 and FR2 on 
large scales (Strom et al. 1981, Barthel et al. |1985 ) with a 
component coincident with the nucleus of the galaxy and 
two lobes of quite different morphology, one (SE) narrow 
and edge-brightened (FR2-like) and clearly connected to 
the nucleus by a low surface brightness jet visible in WSRT 



92 cm data (Mack et al. 1997), and the other (NW) more 



diffuse and centre-brightened (FRl-like) and extended to- 
wards the nucleus but not obviously by a narrow jet. Spec- 
tral index st udies at wavelengths from 92 to 2.8 cm by 
Mack et al. (199S) show that the spectral index in the 
NW and SE lobes gradually steepens from the edge of 
each lobe back towards the nucleus; in addition, the spec- 
tral index in the NW lobe is generally steeper than in the 
SE lobe. 

About two-thirds of the total luminosity originates in 
the steep spectrum (a = —0.6, S oc v a ) nuclear compo- 
nent. Earlier high re soluti on studies (Fomalont & Miley 
Schilizzi et al. 



1975 



1979) have shown its overall size to 



be 1.3 arcsec corresponding to 2 .2 kp c, with an asymmet- 
ric morphology (Barthel et al. 1985|) which is similar in 
some respects to that on the larger scale. 

An HI absorption line has also been detected towards 



the nucleus of 3C236 by van Gorkom et al. (1989) with the 



VLA and mapped in detail with the VLBA by Conway & 



Schilizzi (2000). The high angular resolution results show 
that most of the absorption takes place about 400 pc from 
the nucleus at the tip of the SE jet, and is possibly due to 
a jet-cloud interaction. 

HST/WFPC-2 observations of 3C236 by de Koff et al. 
( |2000D revealed the presence of a dust ring in the nucleus 
~ 10 kpc in diameter and oriented perpendicular to the 
overall radio axis. In addition, a dust lane was observed 
on scales of 1 kpc seemingly connected to the dust disk 
by a strand of dust going from one to the other. 3C236 
is unusually bright at FIR wavelengths (Hes et al 



1995 



Willott et al. 2000| ) which implies that there is warm dust 
in the interstellar medium, presumably the dust seen at 
optical wavelengths. Recent paper by O'Dea et al. ( |2001 ) 
presents HST STIS/MAMA near-UV observations of the 
nucleus of 3C236. Four star- forming regions are found in 
an arc along the edge of the dust lane. 



The age of 3C236 is quite uncertain. An analysis of 
the energy and magnetic field content of the large scale 
lobes by Mack et al. ( [1998 ) led to a maximum age for 
particles in the source of 2 x 10 7 years (for the hotspot 
at the end of the SE lobe). On the other hand, assuming 
an average propagation speed for the hotspot in the SE 



lobe of 0.03c (Lacey & Cole |1993| Scheuer |1995| ), the age 
of 3C236 is 2 x 10 s years. Reacceleration at the hotspot or 
resupply of electrons to the hotspot via a jet is therefore 
necessary unless the speed of propagation is an order of 
magnitude greater than 0.03c. The age of the CSS nuclear 
component is equally uncertain, from a few 10 3 years to a 
few 10 years old, depending on the speed of propagation. 
The estimated ages of the star-forming regions found by 
O'Dea et al. fl200l| ) range from - 10 7 to 10 s - 10 9 years. 

In this paper we report on an extensive study of the 
nuclear component in 3C236 using combinations of VLBI 
and MERLIN data at 18 and 6 cm. Spectral index distri- 
butions of the various components on angular scales of 10 
mas are presented and the physical characteristics of the 
nuclear emission are discussed. We re-examine the HST 



data by de Koff et al. (2000) and its relation to the radio 
structure on both the large and small scales. 



2. Observations and data reduction 

2.1. VLBI data 

Global VLBI observations of 3C236 were carried out 
at 1.663 GHz in September 1984 using telescopes at 
16 locations: Jodrell Bank (25 m Mk2 telescope), Cam- 
bridge (15 m), Defford (25 m), Westerbork-1 (one 25 m 
telescope from the array), Effelsberg (100 m), On- 
sala (25 m), Simeiz (22 m), Iowa (18 m), Haystack (36 m), 
Maryland Point (25 m), Green Bank (43 m), Pentic- 
ton (26 m), VLA-27 (twenty seven 25 m telescopes), Fort 
Davis (25 m), Owens Valley (40 m), and Hat Creek (25 m). 

Additional global observations at 4.987 GHz were car- 
ried out in September 1989 using telescopes at 9 lo- 
cations: Onsala (25 m), Westerbork-1 4 (fourteen 25 m 
telescopes), Jodrell Bank (25 m Mk2 telescope), Effels- 
berg (100 m), Haystack (36 m), Green Bank (43 m), VLA- 
27, Pie Town (25 m), and Owens Valley (40 m). Data 
for the 1.663 GHz observations were recorded with Mk2 
equipment (1.8 MHz bandwidth) and correlated on the 
CIT/JPL data processor at the California Institute of 
Technology. The 4.987 GHz data were recorded with Mk3 
equipment in mode B (28 MHz bandwidth) and correlated 
at the Max-Planck-Institut fur Radioastronomie in Bonn. 



2.2. MERLIN data 

An important addition to the aperture plane coverage for 
3C236 was made by adding MERLIN data obtained at 
the same frequency but at a different epoch to the VLBI 
data. Observations at 1.658 GHz were made with seven 



R.T. Schilizzi et al.: VLBI, MERLIN and HST observations of 3C236 






200 -300 



0.2 0.0 -0.2 -0.4 -0.6 -0.8 -1.0 

ARC SEC 



Fig. 1. The kilo-parsec scale central component of 3C236 
at 1.663 GHz; restoring beam is a circular Gaussian 25 
mas FWHM; contour levels are -1, 1, 1.5, 2.5, 3.5, 5, 10, 
15, 20, 30, 40, 50, 60, 75, 90% of the peak brightness, 
which is 251 mJy/beam. Component B2 is the nucleus. 



Fig. 2. Components B1-B2-C: the central 650 pc of the 
nucleus of 3C236 at 1.663 GHz. The restoring beam is a 
circular Guassian 10 mas FWHM; contour levels are —1, 
1, 1.5, 2.5, 3.5, 5, 10, 15, 20, 30, 40, 50, 60, 75, 90% of the 
peak brightness, which is 216 mJy/beam. Component B2 
is the nucleus. 



MERLIN telescopes in September 1993, and at 4.996 GHz 
with six telescopes in September 1991. In combining the 
VLBI and MERLIN datasets, we assume that the struc- 
ture does not vary between the epochs of observation, and 
that the small scaling errors caused by combination of data 
at slightly different frequencies is acceptable. In particu- 
lar, we assume that the 1-2% variations in the flux density 
of the central compact component seen at 5 GHz over a 
period of 8 years (see sect. 3.2) have a negligible effect 
on the matching of the MERLIN and VLBI data taken 
2 years apart. The phase centres of the VLBI and MER- 



LIN dataset were originally matched (Schilizzi et al. 1988) 



by assuming that the peak emission in the MERLIN im- 
age coincided with component B2 in the VLBI data (see 
sect. 3). There were no common baselines to verify this 
assumption. Howeve r, the VLB A image at 1292 MHz of 
Conway & Schilizzi ( 2000 ) showed that this assumption is 
incorrect and that the phase centre of the MERLIN data 
should coincide with the centre of component C in the 
VLBI data. This has been rectified in the data reduction 
reported here. 

The data were reduced using standard procedures in 
the AIPS package. Images of varying resolution were made 
by applying tapers to the data in the uu-plane. 



3. Results 



3.1. The structure of 3C236 at 1.663 GHz 

Fig. [l] shows the 2 kpc radio structure in the nucleus of the 
galaxy with an angular resolution of 25 mas, and Fig. ^ the 
structure in the central 650 pc with an angular resolution 
of 10 mas. The visibility data derived from the images 
shown are a good fit, within the errors, to the observed 
data on all baselines. 

Fig. [l] shows two main emission areas, the central 
strong emission region (B1-B2-C) and the extended ra- 
dio lobe in the northwest (A) which appears to be a dif- 
fuse plateau structure with low surface brightness. Com- 
ponents Bl, B2, and C are more clearly defined in Fig. |[ 
The brightest and most compact component, B2, was ear- 
lier identified as the location of the nucleus of 3C236 on 
the basis of its spectral index and compactness (Schilizzi 
et al. 1988). Component Bl is a slightly resolved knot in 
the jet to the NW, and component C appears to be a lobe 
with a bright spine which defines the radio jet to the SE. 
The overall radio axis defining the long term orientation of 
the source is given by the line joining the prominent fea- 
tures in the large scale extended lobes and passing through 
the nuclear component; this has a position angle (P-A-) of 
123° ± 2° (Schoenmakers et al. |2000| see also Fig. |). Not 
all the individual components in the 2 kpc radio struc- 
ture lie on this line. One that does is the slightly resolved 
component Bl which lies 77 mas north-west of B2 in P.A. 
— 56?5±0?5. In contrast, the main ridge line of the lobe- like 
component C on the other side of B2, lies in P.A. 116°±2°. 
The most prominent feature in component A is about 720 
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Fig. 3. A slice along the ridge of jet emission seen either side of component B2 at 1.663 GHz in Fig. 2. The profile 
shows that the jet brightness is very similar on either side of B2 as expected for a source near the plane of the sky 



mas from B2 and located north of the overall source axis. 
Weak emission features are seen between this feature and 
the B1-B2-C complex which probably delineate the path 
of the jet to the north-west. Closer to the nucleus, B2, 
the jets appear to oscillate in direction; the initial direc- 
tion of the jets emanating from B2 to the northwest and 
south-east is in P. A. 114° ± 0?5 before changing smoothly 
to a P.A. of 122° about 50 mas from B2 and back again 
to 114° mid- way through component C. These observa- 
tions suggest that the jets emanating from B2 oscillate 
in direction within a cone of half angle ~ 5° — 7°. Note 
that the brightness of the jets close to B2 are very similar 
(Fig. [3]). The transverse structure in component C which 
lies 137 mas south east of B2 may be evidence of a shock 
propagating back through the lobe emission (Barthel et 
al. |1985| , Schilizzi et al. [1988^ Conway and Schilizzi |OOCj) . 



3.2. The structure of 3C236 at 4.987 GHz 

Fig. ^ shows the central radio structure of 3C236 at 4.987 
GHz at an angular resolution of 25 mas, and Fig. || the cen- 
tral 650 pc at a resolution of 10 mas. The main features at 
4.987 GHz match those at 1.663 GHz including the weak 
feature in the north of component A. The surface bright- 
ness sensitivity at 4.987 GHz is not sufficient to recover 
the low level emission in the outer regions of component 
C seen in Fig.s 1 and 2, or to reliably recover the details of 
the low level emission in component A. The double-sided 
jet structure associated with B2 at 1.663 GHz is appar- 
ently one-sided at 4.987 GHz pointing northwest in the 
same direction as Fig ^[ The oscillating behaviour of the 
jet in component C is clearly visible within a cone of half 
angle 5° centred on B2. 

When the full angular resolution of the VLBI data 
is used (1.5 mas, uniform weighting), a different picture 
emerges (Fig. ||). The centre of component B2 is revealed 
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Fig. 4. The kilo-parsec scale central component of 3C236 
at 5 GHz. The restoring beam is a circular Gaussian 25 
mas FWHM; contour levels are -1, 1, 1.5, 2.5, 3.5, 5, 10, 
15, 20, 30, 40, 50, 60, 75, 90% of peak brightness, which 
is 272 mJy/beam. 



as a multi-component structure in position angle 112°, 
with the brightest component at the eastern end. The 
presence of component Bl is also visible in the data al- 
though it is well resolved. This is the highest angular res- 
olution data on the nucleus itself, and it is remarkable that 
at this resolution, the structure appears to be a one-sided 
jet whereas, on larger scales, the structure is two-sided. 

Comparison of visibility amplitudes of these 5 GHz 
data with the corres pondi ng data on 5 baselines from April 
1981 (Barthel et al. 1985 ), shows evidence for brightening 
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Fig. 5. Components B1-B2-C; the central 650 pc of the 
nucleus of 3C236 at 5 GHz; the restoring beam is a circular 
Gaussian 10 mas FWHM; contour levels are —1, 1, 1.5, 2.5, 
3.5, 5, 10, 15, 20, 30, 40, 50, 60, 75, 90% of peak brightness, 
which is 232 mJy/beam. 



of the most compact component in the structure by a few 
tens of mJy over the 8 year period. Comparison of the 
apparent positions of the minima in the visibility data 
shows no consistent evidence for change in the positions of 
individual components at a level of 0.065 mas/year which 
corresponds to an upper limit on the speed of separation 
of 0.4c. 
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Fig. 6. The central 150 pc of 3C236 at 5 GHz; component 
Bl is to the northwest and component B2 (containing the 
nucleus) is to the southeast. The restoring beam is 1.5 
mas; contour levels are 1, 2, 3, 8, 12, 20, 40, 60, 90% of 
peak brightness, which is 85.54 mJy/beam. 



3.3. Source component parameters 

Table 1 gives information on the source parameters mea- 
sured at 1.663 and 4.987 GHz. 1) observing frequency; 
2) total flux density; 3) component name; 4) flux den- 
sity of the component; 5) distance from B2; 6) P.A. of 
the component with respect to B2; 7) mean deconvolved 
component sizes at half power, wl and w2, determined 
from two-dimensional Gaussian fits for the components 
that could be approximated well by a Gaussian; 8) P.A. 
of the component in degrees. The flux densities in columns 
2 and 4 were measured from the images using the AIPS 
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task TVSTAT, and have typical errors of 5%. The sepa- 
rations and position angles in columns 5 and 6 were mea- 
sured from the images by hand and have typical errors of 
5% and 1° respectively. The sizes and P.A.s of individual 
components in columns 7 and 8 were determined by 2- 
dimcnsional gaussian fits where possible and have typical 
errors of 10% and 5° respectively. 

3.4- The spectral index distribution between 1.663 and 
4.987 GHz 

We have constructed spectral index distributions for 
3C236 at angular resolutions of 10 and 25 mas on the 
assumption that the source structure did not vary during 
the time spanned by the observations reported here. There 
is no evidence for flux density variations greater than a few 
tens of mJy (see sect. 3.2) at either frequency. The images 
at 1.663 and 4.987 GHz were convolved to the same res- 
olution, then super-imposed assuming the position of the 
peak brightness in component B2 is the same at each fre- 
quency. Spectral indices were computed for each pixel for 
which the intensities in the two images were both greater 
than 8a. 
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Fig. 8. The spectral index distribution across the central 
650 pc of 3C236 at an angular resolution of 10 mas, for the 
frequency range from 4.987 to 1.663 GHz. The greyscale 
range extends from —1.5 (faintest grey) to 0.5 (black); con- 
tours are shown for spectral indices of —0.9, —0.5, —0.1, 
0, 0.2 
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Fig. 7. The spectral index distribution across the central 
components of 3C236, at an angular resolution of 25 mas, 
for the frequency range from 4.987 to 1.663 GHz. The 
greyscale range extends from —2.0 (faintest grey) to 0.9 
(black); contours are shown for spectral indices of —0.85, 
-0.5, 0, 0.2 

Fig. |^ shows the overall spectral index distribution 
across the source at an angular resolution of 25 mas. The 
spectrum is steep (—1.5 < a < —0.5) over most of the 
structure in the central 2 kpc with the exception of com- 
ponent B2 which has a flat spectral index (a = 0.05) and 
the warmspot in the north of component A which appears 
to have an inverted spectrum (a = 0.3). The other re- 



gions of apparent inverted spectrum in component A may 
not be reliable in view of the substantial difference in uv- 
coverage at the two frequencies; these differences are ca- 
pable of causing ambiguities in the reconstruction of the 
total intensity of this very extended component. It is in- 
teresting to note that the oscillating jet to the SE appears 
to have a flatter spectral index where it passes through 
the centre of component C. 

A more detailed picture of the distribution of spec- 
tral index in the central 500 pc is shown in Fig. g, The 
spectral index in B2 varies from —0.5 in the east to 0.4 
at the base of jet to the northwest; further out along the 
jet it steepens again reaching —1.0. There is also a gra- 
dation of spectral index in component Bl from east to 
west, changing from —1.0 to 0. As remarked earlier, the 
ridge in component C has a spectral index of about —0.5 
in the centre, which falls off to —1.0 towards the edges. 
At the location of the transverse feature in component C, 
the spectral index appears to flatten from —0.9 to —0.5. 
At the eastern end of C, the spectral index of the peak 
seen in the total intensity maps has a flat spectral index 
of about zero. From this point to the leadin g edge of the 
jet is the region in which Conway & Schilizzi ( 200C| ) detect 
maximum HI opacity. 

4. Discussion 

4-1. A re- examination of the HST data 

The HST snapshot of 3C236 made with the WFPC-2 us- 
ing a red broad-band filter centred near 7000 A, shows 
clear evidence of a substantial dust ring (de Koff et al. 
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200C). Fig. H (middle) reproduces the absorption model 
for the galaxy derived by de Koff et al., in which the main 
component of dust is in the form of a ring of radius ~ 
5 kpc whose apparent symmetry axis differs in position 
angle by 15 - 20 ° from the overall radio axis defined by 
the outer edges of the large scale structure (Fig. ||, bot- 
tom). The symmetry axis of the ring is aligned with the 
apparent minor axis of the galaxy itself, within the errors. 
A dust feature ~1.5 kpc long lies internal to the ring and 
approximately parallel to it; this may be the remnant of 
another ring. A further faint dust feature can also be seen 
emanating from the nuclear region towards the SE. 

From the ellipticity of the disk, we derive an apparent 
inclination angle of the radio source to the line of sight of 
~ 60° considering that the radio jets are approximately 
perpendicular to the dust ring. Kotanyi & Ekers (1979) 
first pointed out that large scale radio jets are roughly 
perpendicular to dust disks and lanes, and this has been 
amply confirmed and extended to the nuclear regions by 
de Koff et al. ( |2000D , in particular for well defined disks 
and lanes. It is interesting to note that in NGC 4261, Jaffc 
et al. fllggg ) found a very similar offset in P.A. between 
the radio jet and the symmetry axis of the dust disk in 



that galaxy. These authors refer to Rees (1984) in pointing 



out that Lense-Thirring precession will force the part of 
the accretion disk nearest to the black hole to be axi- 
symmetric about the rotation axis of the hole, whatever 
its spin axis at large radius. At small distances from the 
black hole the jet and disk axes should coincide. In the 
case of 3C236 the linearity of the large scale radio jet to 
the south-east implies that the black hole axis has been 
oriented at 122 ° for the lifetime of the source, and that 
the dust ring axis is offset from this on kpc scales. This 
offset could arise if the dust has been captured from a 
smaller galaxy which has been cannibalised by 3C236. The 
oscillation in direction of the compact and large scale jets 
(see next section), if due to oscillation in the black hole 
axis, may well be related to a merger event. The mass 
of dust in the ring estimated by de Koff et al. (2000) of 
3 x 10 6 Mq seems unlikely by itself to cause the oscillations 
through torque on the accretion disk. 

The inclination of the radio axis in 3C236 implies that 
the northwest jet is approaching and the SE jet receding. 
It also implies that the true size of 3C236 is 4.5 Mpc rather 
than 3.9 Mpc calculated for an inclination angle of 90°. 

The nominal centre of the galaxy in the optical absorp- 
tion image lies about 0.2 arcsec south of the northwest 
end of the small dust feature which extends ~ 500 pc in 
P.A. 110° — 115° before apparently connecting into an- 
other dust strand at P.A. ~145 °. Since the astrometric 
accuracy of the HST images is ~ 0.1 arcsec, it is possible 
that the position of the nucleus coincides with the north- 
west end of this dust feature. Assuming this to be the case, 
the fact that the compact radio jet/lobe to the SE (Fig. 
H top; see also Fig. [j], component C) has the same orienta- 
tion (P.A. ~ 116°) and almost the same extent (~ 370 pc) 



as the feature, suggests that the radio emission and dust 
may be related. 

We can make some qualitative remarks on this ap- 
parent connection. It seems unlikely that the jet has en- 
trained dust on its way out for the simple reason that 
the jet itself is not in contact with the neutral interstellar 
medium. It propagates through the relativistic plasma in 
the centre of the radio 'lobe' which is the shocked plasma 
from the jet that has passed through the hotspot. Indeed 
De Young ( 1998| ) has pointed out that the passage of a 
strong high-speed shock associated with a radio jet will 
destroy dust grains by sputtering processes on timescales 
of 10 4 — 10 5 years, comparable with the age of component 
C. The most the radio jet can do is concentrate the dust 
and hence increase the column density and obscuration 
along certain lines of sight. This could take place via the 
expansion of the radio lobe that generates a bow-shock 
propagating out into the neutral interstellar medium. Of 
course the speed of the bowshock should not be so high 
that the dust is destroyed. 

Perhaps the simplest explanation of the apparent 
radio-dust connection is that the dust feature represents a 
dense cloud that has happened to wander into the path of 
the jet as shown by the HI absorption data. In this case, 
the optical dust feature traces the mass which shapes the 
radio structure, rather than the radio structure generating 
the optical structure via entrainment. 

4-2. Radio morphology 

The radio morphology of 3C236 shows structure from less 
than a few milli-arcseconds to 39 arcmin, a range of 10 6 . 
At first sight, the small and large scale structures are sim- 
ilar in some respects but not in others. The components to 
the NW on the small and large scale have approximately 
the same opening angle, 20° seen from the nucleus, and 
they are both much more resolved in the transverse direc- 
tion than the components to the south-east. The structure 
to the SE is clearly very well collimated on the large scale, 
but on the small scale the SE component flares about 45 
mas (74 pc) from the nuclear component. The higher res- 
olution VLBI data on this component shows that there 
is still collimated flow through the centre of the flaring 
component. The flaring is most likely related to a back- 
flow/transverse shock resulting from the putative collision 
between the jet and a cloud ~ 225 mas (370 pc) from the 
nucleus (Conway & Schilizzi |200C ) . 

The arm-length ratios (the ratio of the distances from 
the nucleus to the extremities of the components) of the 
inner and outer structure are quite different. The end of 
the NW component on the small scales is considerably 
further from the nucleus than the end of the SE compo- 
nent, whereas on the large angular scales, the reverse is 
true. The arm-length ratio for the large scale structure 
is ~ 1.6, defined as the distance of the SE hotspot from 
the nucleus divided by the distance of the end of the NW 
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structure from the nucleus. This is an average arm-length 
ratio for giant radio galaxies (Schoenmakers et al. 2000| ). 
Defined in the same way, the arm-length ratio for the com- 
pact structure is ~ 0.3. The difference may reflect the fact 
that the large scale structure is shaped by the interaction 
of the jet with the inter-galactic medium and the small 
scale structure by the interaction of the jet with the ISM. 
It is also intriguing that the Mpc SE counter-jet is con- 
siderably longer and better collimated than the NW jet. 
This suggests that the jet to the NW has encountered a 
denser galactic or intergalactic medium which has reduced 
its outward motion. The large scale arm-length asymme- 
try c annot be due to light travel time effects (Longair & 
Riley 1979 ), in particular since we conclude that the NW 
side is approaching. 

The fact that 3C236 is so large and has remained so 
stable in overall orientation argues that it resides in a rel- 
atively empty environment. AGN in more dense environ- 
ments (eg. Cygnus A) radiate more of their energy and/or 
become Wide Angle Tail sources, whereas more of the en- 
ergy in objects like 3C 236 g oes into outward expansion 
(eg. Barthel & Arnoud |l996| ). 

The 92 cm image of 3C236 (Fig. § bottom) provides 
strong evidence that the jet to the SE has been on for 
most of the life of the radio source. At the resolution of 
the measurements (45 arcsec), the jet appears to originate 
in the nucleus and continue for ~ 2.5 Mpc before termi- 
nating in a double hotspot (Barthel et al. |1985| ). The jet 
appears to broaden in the transverse direction as it moves 
further away from the nucleus, eventually subtending an 
angle of ~ 10° which is considerably larger than the angle 
subtended by the s epara tion of the hotspot at the nucleus, 
3°. Barthel et al. ( 1985 ) show that the ridge line of the 
SE structure wiggles over the final 600 kpc of its trajec- 
tory, so that the cone of opening angle 10° may well be 
the envelope of the jet's transverse motion. The jet ridge 
line in the 92 cm image is located along the southern edge 
of the cone but curves slightly northward before reaching 
the double hotspot. On average, the jet and counterjet di- 
rections have remained very stable over the lifetime of the 
source. 

Close to the nucleus, the SE jet is apparently blocked 
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due to a collision with a cloud (Conway & Schilizzi 2000) 



which has caused what looks like a transverse reverse 
shock back towards the nucleus. This is presumably of 
recent origin since the jet must re-start further out, as we 
see in the 92 cm data. The highest resolution da ta on the 
large scale jet was published by Barthel et al. ( 1985 ) at 
21 cm but does not have sufficient sensitivity to trace the 
jet in close to the nucleus. 

Superficially, the evidence from the large-scale jet to 
the SE is that activity in 3C236 is continuous rather than 
recurrent. However, the 92 cm image in Fig. ^| has a res- 
olution of 45 arcsec or 74 kpc projected on the plane of 
the sky. Gaps in the emission of this size would not be 
detected. Taking a jet speed of 0.1c (a typical speed for 



& Conway 1998; Tschager et al. 200C), and assuming that 
we actually see the jet and not the backflow in a thin lobe, 
we would not detect 'off' periods in jet activity lasting less 
than 2 x 10 6 years. This period is longer than the probable 
length of the present 'on' period represented by the kpc 
central component, 10 4 — 10 5 years. The evidence from the 
structure to the NW points more towards episodic rather 
than continuous activity. In the central kiloparsec, the jet 
is traceable via a number of knots on its way from the 
nucleus via component Bl out to the warmspot in compo- 
nent A at a distance of ~ 1.2 kpc. Insufficient sensitivity 
may limit our ability to detect more of the jet on these 
sub-galactic scales. On larger scales, the emission extends 
~ 1.4 Mpc in what appears to be a more relaxed jet than 
its well-collimated counterpart to the SE. In the 21 cm 
image made by Barthel et al. (1985) there is a prominent 
feature in the NW structure which suggests that the en- 
ergy supply may have been variable. At 92 cm (Fig. |^ 
bottom) and at 49 cm (Mack et al. 1997), there is emis- 
sion south and west of the main radio axis to the NW. 
The origin of this emission is unclear. 

Double-double radio galaxies (Kaiser et al. 2000D ap- 
pear to provide good evidence that activity in AGNs can 
be recurrent. Presumably these are objects in which the 
'on-off ' cycles are comparable. There may be a wide range 
of duty cycles amongst AGNs; 3C236 may be an example 
of a source which is occasionally 'off', perhaps due to the 
jet channel collapsing or the flow being blocked by interac- 
tion with a cloud or a drop in supply of fuel to the central 
black hole. The central structure in 3C236 may reflect the 
normal level of activity in the nucleus. 

Mack et al. ( 1998| ) and Schoenmakers et al. ( [2000 ) have 
estimated the age of 3C236 from synchrotron loss models 
and spectral aging as ~ 10 8 years. The latter authors also 
estimate the speed of advance of the western jet into the 
intergalactic medium as 0.1 c. This speed is at the lower 
limit of what can be measured with VLBI for compact jets 
in the nu cleus (eg. Owsianik & Conway 199S, Tschager et 
al. |2000|) . 



Table 2 lists physical parameters for the core compo- 
nents of radio source 3C236 calculated assuming equipar- 
tition, ellipsoidal geometry for components Bl, B2 and 
the NW lobe, cylindrical geometry for the component C, 
and unity filling factor. The radio window has been taken 
from 10 7 — 10 11 Hz. Component sizes have been measured 
from the 5 GHz images. The volumes derived for compo- 
nents A and C may be in error by a factor of 2-3; for Bl 
and B2 the errors are of the order of 10%. The spectral 
indices in the Table are mean values for the components, 
with typical errors of 10%. The radio luminosity and en- 
ergy in particles and fields have been calculated following 



Pacholczyk (1970), assuming the proton-electron energy 
ratio k — 10. The estimated errors in these quantities 
vary with the component in question: for components A 
and Bl the values listed may be in error by a factor of 2 
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to 3, while for components B2 and C the error estimates 
are more like 10 to 30 %. The parameter values found are 



comparable with those deduced by Barthel et al. (1985) 
taking into account differences in the initial assumptions 
and measured component parameters. 

The jet power measured on large scales by Schoen- 
makers et al. (2000) is ~ 10 45 ergs -1 assuming an upper 
limit of 100 GHz to the radio window. In contrast, on the 
smaller linear scales, the jet power is ~ 5 x 10 43 ergs -1 
with the same upper limit to the radio emission. This may 
indicate that the nucleus was more active in the past when 
the outer lobes were created. 

Can we reconcile the symmetric structure seen on 
10 — 20 mas scales at 1.663 GHz with the one-sided struc- 
ture seen at 5 GHz on 1 mas scales? The HST image (Fig. || 
middle) shows convincingly that 3C236 is ~ 30° from the 
plane of the sky with the NW jet approaching. For radio 
emitting plasma moving along the jets with a bulk mo- 
tion less than a few tens of percent of the speed of light, 
no beaming effects should be detectable. This is consistent 
with the symmetric morphology seen in the nucleus of the 
galaxy on 10 - 20 mas scales (see Fig. ||). The direction of 
the one-sided structure on milli-arcsecond scales is consis- 
tent with the NW jet being the approaching one, but the 
brightness of the first component in the jet is ~ 50 times 
stronger than the noise in the receding jet direction. This 
is not consistent with the angle to the line of sight being 
60°, independent of the magnitude of the bulk velocity 
(Baum et al. 1997). A boost in intensity of a factor of 50 
can be achieved if the compact jet is oriented at 35° to 
the line of sight rather than 60°, and the bulk velocity 
is > 0.8c. This implies a substantial change in the orien- 
tation of the small scale jet which is perhaps unlikely in 
view of the long term stability of the large scale jet sys- 
tem. However, the P.A. of the compact jet on the plane 
of the sky is 10° different to the large scale P.A., indicat- 
ing that the jet orientation does change. A more thorough 
search for superluminal motion in this structure may be 
illuminating. 

An alternative possibility is that the lack of a coun- 
terjet on parsec scales is caused by free-free absorption by 
ionised gas in the inner region of the gas-dust disk sur- 
rounding the black hole. This is seen on VLBI scales in 
NGC 4261 (Jones et al. p000|) Cen A (Jones et alpgg ), 
and 3C84 (Walker et al. |1994| ; Vermeulen et al |1994| ). In 
NGC 4261, the 1.6 GHz image is symmetric about the nu- 
cleus as is the case for 3C236, but at 8.4 GHz NGC 4261 
shows the effects of free-free absorption. Higher frequency 
observations are required to test this possibility in 3C236. 



5. Conclusions 

The radio structure in the nucleus of the giant radio 
galaxy, 3C236, has been investigated using combinations 
of VLBI and MERLIN data at 1.663 and 4.987 GHz. The 
nuclear structure is dominated by two main complexes of 



emission whose orientation and separation are in general 
well aligned with the large scale structure, although differ- 
ences in the projected P.A. reach ~ 10° within a few tens 
of milli-arcseconds either side of the nuclear component 
indicating that the jets oscillate about the mean direction 
set by the Mpc scale structure. The nuclear component 
has been identified by its flatter spectrum and compact 
size in spectral index distributions made at angular reso- 
lutions of 10 and 25 mas. The emission is one-sided to the 
north-west on scales of a few milli-arcsconds, but symmet- 
ric on scales of a few tens of milli-arcseconds. 

Re-exa minati on of HST data on 3C236 taken by de 
Koff et al. ( [2000 ) reveals that the normal to the plane of 
the dust ring in the centre of the galaxy is at ~ 60° to the 
line of sight and parallel, within 15° to 20°, to the pro- 
jected elongation of the radio structure. Assuming that 
the radio axis is parallel to the normal to the dust disk, 
we identify the shorter north-west jet as being on the ap- 
proaching side, and the longer south-east jet as receding. 
The true size of 3C236 is therefore 4.5 Mpc rather than the 
3.9 Mpc determined assuming the source is in the plane 
of the sky. Coincidence, within the astrometric errors, of a 
dust feature with the compact south-east jet (component 
C) is unlikely to be evidence that the south-east radio jet 
entrains dust in the first 400 pc after leaving the nucleus, 
rather that the jet has encountered a cloud of gas and dust 
which actually shapes the radio component. 

The one-sided nuclear radio structure to the north- 
west can be explained on the basis of relativistic beaming 
effects, if the angle to the line of sight is ~ 15° and the 
bulk motion of the radio emitting plasma is highly rela- 
tivistic. However, it is not clear that such a large deviation 
in angle to the line of sight compared to the large scale 
structure can be reconciled with the relatively small de- 
viations in projected orientation of the jet system from 
small to large scales. If the structure is due to beaming, 
then we must conclude that the change from one-sided 
to two-sided emission implies that the angle to the line of 
sight has reached 60° at a separation of ~ 20 mas from the 
nucleus. The alternative possibility of free- free absorption 
of the counter jet needs investigation. 

At the resolution of WSRT data at 327 MHz by 



Schoenmakers et al. ( 2000 ), the south-east jet appears es- 
sentially continuous over a distance of 2.5 Mpc from the 
nucleus to the hotspots, making this the largest jet known 
in the universe. It is likely, however, that the central ac- 
tivity is episodic but with a shorter duty cycle than for 
the double-double galaxies studied by Schoenmakers et al. 
The current central radio emission in 3C236 is evidence of 
a recent major burst of radio emission superimposed on 
a more uniform level of activity, perhaps as a result of 
subsuming a companion galaxy. 
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Table 2. Component parameters of 3C236 



Components 


Volume (cm' 1 ) 


Sp index 


Lradio (ergs x ) 


E e (erg) 


E H (erg) 


A 


10 b2 


-0.4 


5 x 10 4i 


1Q b3 


8 x 10 b4 


B2 


< 10 59 


+0.05 


5 x 10 42 


1Q 53 


8 x 10 53 


Bl 


10 59 


-0.6 


> 2 x 10 41 


4 x 10 52 


> 3 x 10 53 


C 


10 61 


-0.9 


4 x 10 41 


7 x 10 53 


6 x 10 54 
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Fig. 9. The images of 3C236 from observations with VLBI 
at 18 cm (top, see Fig. 1 for details), HST absorption 
at 7000 A (middle) and the WSRT at 92 cm (bottom, 
Schoenmakers et al. |200C ). For the WSRT map, the restor- 
ing beam is 45 arcsec; contour levels —0.1, 0.1, 0.3, 0.5, 
0.8, 1.4, 2.4, 8, 15, 30, 40, 60, 90% of peak brightness, 
which is 7.24 Jy/beam. The optical nucleus is at (0,0) in 
the middle panel. 



